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Number-conserving (or conservative) cellular automata have been used in several 
contexts, in particular traffic models, where it is natural to think about them as 
systems of interacting particles. In this article we consider several issues concern- 
ing one-dimensional cellular automata which are conservative, monotone (specially 
"non-increasing"), or that allow a weaker kind of conservative dynamics. We intro- 
duce a formalism of "particle automata", and discuss several properties that they 
may exhibit, some of which, like anticipation and momentum preservation, happen 
to be intrinsic to the conservative CA they represent. For monotone CA we give 
a characterization, and then show that they too are equivalent to the correspond- 
ing class of particle automata. Finally, we show how to determine, for a given CA 
and a given integer b, whether its states admit a 6-neighborhood-dependent rela- 
belling whose sum is conserved by the CA iteration; this can be used to uncover 
conservative principles and particle-like behavior underlying the dynamics of some 
CA. 
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1 Introduction 



Cellular automata (CA) are discrete dynamical systems, where states taken 
from a finite set of possible values are assigned to each cell of some regular 
lattice; at each time step, the state of a cell is updated through a function 
whose inputs are the states of the cell and its neighbors at the previous time 
step. They are useful models for systems of many identical elements when 
the dynamics depends only on local interactions. Conservative (or "number- 
conserving") cellular automata represent a special class of CA, in which the 
sum of all the states, that are integers, remains constant as the system is 
iterated. This property arises naturally when modeling phenomena such as 
traffic flow ([NS92]), eutectic alloys ([Koh89,Koh91], or the exchange of goods 
between neighboring individuals. When number-conservation is not apparent 
for the initial system, its detection can be interesting by itself, and may help 
to prove dynamical properties. 

Necessary and sufficient conditions for a CA to be number-conserving are given 
in [BF98] for one dimension and states {0, 1}, and in [BF01] for one dimension 
and states {0,... , q — 1}; a generalization for two and more dimensions is 
found in [DFR03]. In [Mor03] the definition — and the characterizations — are 
extended to allow general sets of states S C Z, and an algorithm is given to 
decide, for any CA, whether its states can be relabeled with integer values, 
so as to make it number-conserving. In [MI98] and [MTI99] the universality 
of reversible, number-conserving "partitioned" CA is proved for one and two 
dimensions, respectively. In [Mor03] the universality of usual (not partitioned) 
conservative CA in one dimension is proved. In fact, it is shown that any 
one-dimensional CA can be simulated by a conservative CA; this proves the 
existence of intrinsically universal conservative CA in the sense defined in 
[O1101]; this notion of universality is stronger than the usual one (the ability 
to simulate universal Turing machines). A construction of a logically universal 
conservative CA in two dimensions is given in [IFIM02] ; they also construct a 
self-reproducing model in a two-dimensional conservative CA, by embedding in 
it the well known Langton's loops. Another interesting work is found in [DFR], 
where the CA classifications of Kurka [Kur97] and Braga [BCFM93,BCFV95] 
are intersected and the existence of conservative CA in the resulting classes 
is checked. A recent article by Fuks [Fuk] considers probabilistic conservative 
CA. 

In the articles of Boccara and Fuks ([BF98,BF01]) the necessary and sufficient 
condition was used to list and study all the conservative CA rules with small 
neighborhoods and small number of states. For all the rules they study, they 
give a motion representation: the state of a cell is interpreted as the number 
of particles in it, and the CA rule is interpreted as an operator that governs 
the interaction of these identical, indestructible particles. Fuks [FukOO] and 
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Pivato [Piv02] have independently shown that this interpretation is always 
possible (in the one-dimensional case). In the same spirit but with a very 
general definition of "particles", Kurka [Kur03] has recently considered CA 
with vanishing particles. 

For the sake of completeness and to avoid confusions, it is worth mentioning 
other contexts in which particles have been considered. On one hand, there 
are the interacting particle systems (IPS), with a long history in probability 
theory [Lig85], and the lattice gases, some of them with associated CA models 
[Boo91]; in general, they cannot be written as conservative CA. The well- 
known two-dimensional Margolus CA [Mar84] is number- conserving and was 
designed to allow rich interactions of particles; it does not fit in the definition 
given here, because of its alternating neighborhood. Particles have been widely 
used in computer graphics [HE88] , sometimes using CA with the neighborhood 
of Margolus [TET95]. The word "particle" is also used to describe emergent 
particle-like structures that propagate in CA [BNR91,DMC94,HC97,HSC01]; 
in this last sense, it is close to the spirit of our last section, and to that in 
[Kur03]. 

In this article we consider one-dimensional cellular automata; Section 2 gives 
the necessary definitions and reviews (and generalizes) some relevant previous 
results, while Section 3 gives our definition of particle automata (PA) as a for- 
malism for motion representation. Section 4 deals with several issues related to 
conservative CA. First we prove (again) their equivalence with the (conserva- 
tive) PA; then we discuss several behaviors that PA may exhibit, showing that 
some of them (like anticipation and global cycles) may be intrinsic to some 
conservative CA. We also consider the special properties of state-conservation 
(where a sensible particle representation will recognize each state as a different 
kind of particle) and momentum preservation (which, in spite of being defined 
in terms of the PA, depends only on the conservative CA it represents). The 
main result of the paper is in Section 5, where we characterize non-increasing 
CA and show how to represent them with particle automata. Finally, Section 
6 considers CA where the states can be relabelled, in a way that depends 
on the neighbors of a cell, in order to obtain a conservative dynamics (and a 
representation in terms of particles). 



2 Definitions and Some Previous Results 

Cellular automata: A one- dimensional cellular automaton (CA) with set 
of states Q = {0, ... , q — 1}, is any continuous function F : Q z — > Q z which 
commutes with the shift. It is well known that cellular automata correspond 
to the functions F that can be expressed in terms of a local function: F(c)i = 
/(ci+jv), for all c G Q z , i G Z, and N a fixed finite subset of Z, called the 
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neighborhood of F. N can always be assumed to be an interval of integers which 
includes the origin, and we write F(c)i + d = /(cj+jv), with iV = {0, ... , n — 1} 
and d G Z an offset; rules with the same / but different d will be identical up 
to a shift. It will be useful to define, for n G N and Q — {0, ... , q — 1}, 

CA(q,n) = {.[:(.)<" " > ■ (.)} 

Any CA can then be expressed by an element of CA(q, n) for some q and n, 
combined with an offset d which tells where the image of the neighborhood is 
placed. The 256 elementary CA, for instance, correspond to £4.(2,3), usually 
with d — 1. £4 will denote the union of CA(q, n) over all q and n. 

A common shorthand notation for cellular automata is the codification used 
by Wolfram [W0I86]: the code for an element / eCA(q, n) is given by 

Code(/) = £ f(xi,...,x n )qELi<r-"*>> 
(xi,... ,x n )eQ n 

Configurations: An element in Q z is called a configuration. A configura- 
tion is said to be finite if all but a finite number of its components are 0. A 
configuration c is said to be periodic if q = Ci +P , for all i, for some p G Z, 
p 7^ 0; in this case, p is said to be a period of c. 

Monotone and conserved quantities: Consider a CA F on Z, and let 

Cp be the set of all periodic configurations in Z; for each c E Cp choose a 
period p(c). Let be a function : — > R, where 6 is a nonnegative integer. 
Then is said to be a non-increasing additive quantity under F if and only if 

p(c)-i p(c)-i 
£ fiiFi^k, - , -F(c)jfc +6 _i) < £ 0(c fc , ... ,c fc+6 _i), Vc G C P . (1) 
fc=0 fc=0 

Similarly, is said to be non- decreasing additive quantity if condition (1) 
holds with the inequality in the other direction. It is easy to see that is 
non- increasing if and only if —0 is non-decreasing. If is both non-increasing 
and non- decreasing, it is said to be a conserved additive quantity (in this case, 
(1) holds with an equality sign). We say that is monotone if it is either non- 
decreasing or non- increasing. In [HT91] it is said that "an additive conserved 
quantity is a discrete-time analog of what we usually call a conserved quantity, 
such as energy, momentum and charge of a physical system" ; the sentence can 
be rephrased for the monotone case. 

Finitary characterization: The previous definitions consider the addition 
of a density function over a period of a periodic configuration. Another pos- 
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sibility would be to consider the sum over finite configurations: we may say 
that is a finitely non- increasing additive quantity if condition (1) holds for all 
finite c in Z, instead of Cp, with the sums being taken now over the whole Z. 
(Here we are assuming that 0(0,... ,0) = 0; if this is not the case, we consider 
= — 0(0, ... , 0) instead.) It turns out that the two notions are equivalent: 

Theorem 1 (Generalized from [DFR03]) Let F be a CA andcj) be a func- 
tion : S b — > ffiL Then is an additive conserved (non-increasing, non- 
decreasing) quantity for F if and only if it is an additive finitely conserved 
(non-increasing, non- decreasing) quantity for F . 

Sketch of the proof. In [DFR03] the equivalence is proved for conserved 
quantities, in dimension 1, when : Q — > Q is the identity; however, their 
proof includes both the non-increasing and the non-decreasing cases, and can 
be easily extended to the case of a general 0. In one direction the proof is 
trivial: if the condition holds for all periodic configurations, and c is a finite 
configuration, then the condition is shown to hold for c by building a periodic 
configuration with blocks that include the non-zero part of c. On the other 
hand, if the condition is not verified by a periodic configuration with repeated 
word w, then it will be not verified for a finite configuration of the form 
... OOOw^OOO... , for N large enough: the surplus (or deficit) of the periodic 
configuration is amplified by the growing N, while the only terms that could 
reduce it (those corresponding to a neighborhood of Ow and wO) remain fixed. 

Notice that the same argument can be also extended to higher dimensions: by 
repeating an n-dimensional pattern enough times, its surplus will be amplified 
as N n , while the terms corresponding to the border, though not fixed, will grow 
only as N 71 ' 1 . □ 

The following theorem is a useful characterization of conserved quantities in 
one-dimensional CA. 

Theorem 2 (Hattori and Takesue [HT91]) Let F be a one- dimensional 
CA with local rule f G CA(q,n). Let a be an arbitrary element in Q = 
{0, ... , q — 1}. Then : Q b — > R is an additive conserved quantity under 
F if and only if 

0/(x O , ... , X b+n ^ 2 ) - (f>(x , ... , X b -i) 



b+n-2 




) + - ,a,Xl, ... , X b+n -i-i)} 



(2) 



b— i b—i 
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for all xq, ... , Xb+n-2 G Q, where 

(f>f(x Q , ... ,x b+n - 2 ) = (f)(f(x ,...,x n - 1 ),...,f(x b - 1 ,...,x b+n - 2 )) 



Monotone and conservative CA: A cellular automaton is said to be non- 
increasing (non- decreasing, conservative) if the identity of its set state is a 
non- increasing (non- decreasing, conservative) quantity for its dynamics. Since 
the condition depends only on the local rule of the CA, and not on the offset, 
we will define CA~(q, n), CA°(q,n) and CA + (q, n) to be the non-increasing, 
conservative and non- decreasing rules in CA(q,n), respectively. 

Theorem 2 implies that / G CA°(q,n) if and only if, for all (x 1: ... ,x n ) G Q n , 

> (3) 



n—l 



f(x l7 ... ,x n ) = x 1 + ' 
k=l 



/((V^O, x 2 , ... , x fc+ i) - /(JV^O, x x , ... , x k ) 

n—k n—k 



This characterization, given in [BF01], can be generalized to higher dimen- 
sions, though its explicit form becomes hard to write (it was done by [DFR03]). 



Some more notation: The letter q will always denote the number of states, 
and the letter Q will denote the set {0, ... , q — 1}. With we will denote an 
infinite sequence of zeroes (thus, OwO denotes a word w surrounded by infinite 
zeroes). For / G CA(q,n), we will denote by f{u/v) the block image of word 
u when followed by word v : 

f( U /v) = f(w , ... ,W n _i)/(wi,... ,W n ) ... ... ,w M+n _ 2 ) 

where w = uv. Of course, only the n — l first elements in v contribute to 
f(u/v). Furthermore, we denote f{u) = f(u/e), for any u with \u\ > n, where 
e is the empty word. 

With this notation, the condition for / G CA(q, n) to belong to CA~(q, n) can 
be restated as 

\w\-l \w\-l 

E f(w/w) k < ]T w k WweQ* (4) 

fc=0 k=0 

If #(•) is a vector-valued function, we will use the notation [#(•)]& to refer to 
its A;-th component. Finally, we will use the function (•)+, defined as (x) + = x 
for x >= 0, and (x) + = for x < 0. 
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3 Particle Automata and Motion Representations 

A common way to look at conservative CA is through their representation 
in terms of particles: the state of each cell is interpreted as the number of 
particles contained in it, and a rule is given describing the motion that these 
particles will have, depending on the local context; we will add the possibility 
of vanishing, and will formalize this as particle automata. 

A particle automaton (PA) with set of states Q = {0, ... , q — 1} will act on 
Q z , like a CA, and like a CA it will be defined by a local rule (or rather, by a 
set of rules) taking as input the states of a local neighborhood, which is again 
some N of the form {—£, ... , r}. 

For a local configuration w = c_^, ... , Co, ... , c r with Co > 0, a function g C() 
will give the new positions of the Co particles at the origin: we have a set 
(9i)i=i,... ,g-i of functions 

9l :Q e xQ r -,(NU]Y 

where the dagger (f) represents the "vanishing" option. Thus, for a configu- 
ration c G Q z ', the k-th particle at position i (with 1 < k < Cj) will 

{vanish if [g Ci (d- t) ... , Cj_i, c i+1 , ... , c i+r )] k = f 
move to position % + [# Ci (q_^, ... , c m , ... , c i+r )\ k otherwise 

Since the particles are undistinguishable, the order of the components of each 
gi is irrelevant; only the number of components mapping to each element of 
(N U f) is important. 

Some examples may help to make the definition clear. Consider a PA with q 
states that moves all particles one position to the left. It will have £ = 1, r = 
(i.e., N = {-1,0}), and 

g a (P) = (-1, ... , -1) for all < a < q, < (5 < q 

S v ' 

a 

As a second example, consider a PA with Q = {0, 1, 2}, I = 1, r = 2, and g±, 
g 2 such that <?i(0, 2,0) = (f) and g 2 (l,0,0) = (0,2). Figure 1 shows its effect 
on the configuration 0120. 

Thus, a PA G is defined by a tuple G = (q, N, (<7i)j=i,... , q -i)- The global action 
of G on c G Q z is given by 

G(c)i = min{ q-1, #{(j, k) : j + [g c ,(cj-e, - , c,_i, c i+1 , ... , c j+r )] fc = «} } 
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Fig. 1. Example of particle motion 

In other words, G(c)i is the number of particles arriving at i, with a maximum 
of q— 1; this last condition prevents "overflows", but, in any case, any rule can 
always be fixed to avoid needing this overflow control, by sending the excess of 
particles to "f" (though this may require an extension of the neighborhood). 

We will denote by VA(q, n) the set of all PA with q states and neighborhood 
of size \N\ = £ + r + 1 = n, and by VA°(q, n) the set of conservative PA: those 
members of VA(q, n) such that their functions gi go to N' 1 (nothing vanishes), 
and avoid the overflow (i.e., the minimum in the above definition is always 
the right term). Notice that unlike the elements of CA(q,n), an element of 
VA(q, n) is not combined with an offset to define the actual PA. 

Motion Representation: In [BF98] and [BF01] Boccara and Fuks give 
motion representations for each of the rules they study; we will follow their 
notation, which is an efficient and intuitive way of expressing a PA. A motion 
representation is a list of specific local configurations of a given a cell, with 
arrows indicating the motion to be performed by the particle(s) located in 
that cell for all these local configurations. As a simplification, if for a given 
local configuration, the particles do not move, then this configuration is not 
listed. Some examples may clarify this point. Consider 

2 1 

M 1 = {10,0011} , M 2 = {20,21} 

Mi, for q = 2, is read as follows. If a particle sees an empty cell on its right, 
it moves to it. If that cell is occupied, but the two cells on its left are empty, 
it moves to the closest one. In any other case, it keeps its current position. 
(The neighborhood in this case is {—2, —1, 0, 1}.) Numbers may be added to 
the arrows when a cell may be occupied by more than one particle. This is the 
case in M 2 , with q — 3. If two particles are in a cell, then one, two, or none 
of them may move to the right, depending on the space available there. Here 
are two more examples: 

k 

M 3 = { k = l,...,q-l} , M 4 = {10} 

In M 3 the bullets (•) are wildcards. In this example, everything moves two 
steps to the left, regardless of what the states in the other cells are. If we 
denote by a the shift (the CA rule such that /(ao, di) = fli, with offset 0), 
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then M3 represents the CA rule a 2 = a o o. Rule M4 (for q = 2) shows that a 
particle will move to the right if, and only if, that cell is empty; its effect on 
{0, 1} Z is the same as the elementary CA rule 184, with offset 1. The notation 
is easily extended to include vanishing particles, by adding a hat (") for each 
of the particles that vanish: 

2 1 

M 5 = {01, 11, 21} 

Here the l's will travel to the left until they meet a 2, and then will disappear. 

Note that PA do not distinguish particles; if k particles arrive at a same 
cell, at the next time step the rule says how many of them will move to 
each neighboring cell, but does not say which particles are moving. However, 
particles are, in a sense, distinguishable, since we know how many go from each 
cell to each other cell at each iteration. To say that "a particle moved from j to 
i, while another moved from i to k" , is not the same as saying "a particle moved 
from j to k, and another stayed at i" . When arbitrarily large groups move, 
as for instance in the shift, we have to assume that some intermediate cells 
with unchanging values are changing the particles they contain; otherwise, we 
would need an infinite neighborhood to describe the motion. If we implement 
the system and want to trace the particles throughout the iterations, we need 
to add some criterion. A sensible choice (implicitly applied in [BF01]) is to 
keep the order of the particles along the line. 

Both CA and PA are mappings from S z into itself. We say that a CA F is a 
projection of a PA G if F(c) = G(c), for all c E S z . The following proposition 
is straightforward. 

Proposition 3 Let G be a PA and F be a projection of G. Then G G VAif 
and only if F e CA' , and G G VA° if and only if F G CA° . 

3. 1 CA for a given PA 

Theorem 4 For any G G VA(q,n), there is a unique CA F with local rule 
f G CA(q, 2n — 1) which is a projection of G. 

Proof. The function G is obviously continuous and shift-commuting. Hence, 
it may be written as a CA, which is uniquely defined but for the neighborhood 
(since the neighborhood size can always be increased). The only thing we have 
to check is that a neighborhood of size 2n — 1 is enough. For this note that, in 
the definition of PA, the state of c[ is completely determined by the particles 
that will move to i. If we denote by {—£, ... , r} the neighborhood of G (n — 
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£ + r + l), these particles have to be in the cells {i — r, ... ,i + £}; their behavior, 
in turn, is completely determined by the values in U*ii_ r {j — £, ... , j + r}, i.e., 
by {c-r-e, ... , c r+e }. □ 

We will denote by 11(G) the smallest CA F (with respect to \N F \) that is a 
projection of G. If, for G, we also take the minimum possible neighborhood, 
then, in most cases, we have \N G \ < \N F \. Take, for instance, the motion rule 
M 4 defined above: rule 184, a 3-input rule, is the smallest CA matching it. 
The reason is that an occupied cell must know if its particle will leave, i.e., 
must look to the left, and an empty cell must know whether a particle will 
arrive, i.e., it must look to the right. The relation \Nq\ < \N F \ is, however, 
not always verified, as shown by the following examples. 

Example 1: Consider the PA G = (2, {—3, ... , 2}, g) with g described by the 
following motion representation: 

M 6 = { 0001 , 1001 , 110 , mo } 

It may be checked that 11(G) = ({0, 1}, {—2, ... , 2}, /), with / given by code 
number 3221127170. A way to look at this situation is this: from the informa- 
tion given by the occupancy numbers of the neighboring cells of an occupied 
cell, we know that the particle will leave it and go to the left by looking at the 
two cells in that direction, but we do not know its precise destination. If we 
take the viewpoint of the destination cells, then we know it, since they do see 
the rest of the (PA) neighborhood of the particle. In other cell knows 

that it will remain occupied, but it does not know if the occupying particle 
will be the same. □ 

Example 2: Here we show a family of PA for which the minimal CA actually 
requires the whole neighborhood allowed by Theorem 4. Consider the family 
of PA Gi )T with q = 3 described by the motion rule 

M 7 = { 2uiv00, 6ulv2 } where u = 2 e ~ 1 and v = 2 r_1 

If we have a local configuration aulvOulvfl, which has size 2£ + 1 + 2r, the 
next state of the cell in the middle, which now contains a 0, depends both on 
a and (3. Hence, the minimal CA has neighborhood {—£ — r, ... , £ + r}, while 
G requires only {—£, ... ,r}. 

If we want an example with only two different states, the PA needs to be a 
bit more complicated; one possibility is described by the motion rule 

M 8 = { OOulvO , ilulvO , lib } where u = e ~ 2 and v = r_1 
□ 
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4 Particles for Conservative CA 

4-1 PA for a given conservative CA 

Theorem 5 For each conservative CA F with local rule f G CA°(q,n), there 
exists G G VA°(q, 2ra - 1), such that 11(G) = F. 



Proof. This result was independently proved both by Fuks [FukOO] and Pi- 
vato [Piv02]. Though very different, all these first proofs (including an unpub- 
lished one by ourselves) are rather long; the idea, however is very simple, and 
we sketch it here. In addition, Theorem 15 in Section 5.4 will generalize it, 
providing yet another proof for the theorem. 

Let {—£,... , t} be the neighborhood of F. We define G with Nq = { — {£ + 
r),... ,(£ + r)}; for all w = w-£- r ,... ,we +r E Q 2n ~ x with wq > 0, we have 
to define the new positions g WQ (w-e- r , ... , W-±, wi, ... ,W£ +r ). This is done by 
computing the image of the configuration c = OwO, with wo at the origin, and 
matching the particles in the image d with the particles in c, from left to right. 

We claim that the new positions of the Wq particles at the origin are in 
{— r, ... , £}. Suppose that one of them goes to a position to the left of — r: this 
means that the particles in (cj)j<o did not match all the particles in (cQj<_ r . 
Since these d i depend only on the values of (cj)j<o, this would contradict / G 
Cv4°for the configuration 0w_^_ r ... u>_i0. The symmetric argument shows that 
no particle from the origin moves to the right of position £. 

Now take any u G Q*, and define c = OuwO and d = F(c). As before, we 
match the particles in c and d, from left to right. Notice that we have added 
I]i=o 1 Ui particles in the preimage; since / G CA° , we have also added the 
same number of particles in the image, and they must be to the left of — r (the 
rest of the image has not changed). Hence, the movement of the wo particles 
at the origin is the same as before, and we see that the effect of applying G 
on any finite configuration c is the same as the effect of matching the particles 
with those from the image of c through F. We conclude that F and G have 
the same effect on c, and F = 11(G). □ 

Corollary 6 In Theorem 5, if the CA F has offset £ (and hence neighborhood 
{—£, ... , r}, with r = n — 1 — £), then G has neighborhood { — (£ + r), ... , (£ + r)}, 
and for each i = 1, ... ,q — 1, gi(Q 2e+2r ) C {— r, ... ,£} 1 (i.e., the particles move 
only to {— r, ... , £} ). 

Example 3: If we take / G O4 (2,5) with code #288146448 (see Table 1) 
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and offset 2, we obtain a PA with a motion representation given by 

M 9 = { U00 , 1001 , 1010 , 1011 , flOl } 
If we set the offset to 0, the result is 

Mio = { «1011 , «1101 , ^11 } 

Example 4: If we take / e £A°(2,4) with code #49024 (see Table 2) and 
use offset 1, the resulting PA has motion representation 

M n = { 1«0 , 100 } 

□ 
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Tabie 1 

Lookup table for rule #2881464448 in CA°{2,5) 
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Table 2 

Lookup table for rule #49024 in £A°(2,4) 
4-2 On Some Behaviors of PA 

The theorem in the preceding section shows how to construct, for a given 
conservative CA F, a conservative PA G such that F = 11(G). We will call 
this the canonical PA for F, since it is the only one that preserves the order 
of the particles along the line. However, it is not the only PA that matches 
F (in fact, there are infinite PA matching any given CA). For this reason, we 
shall discuss some behaviors that a PA may exhibit, and whether or not they 
can be intrinsic to certain CA; this may be relevant in the applications. 

Local cycles: For a PA and a given configuration, we say that there is a 
local m-cycle in the iteration if there is a chain of particles Po?Pi 5 ■•• ,Pm, with 
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Po = Pm, all located in different cells, such that, for i — 0, ... , m — 1, each pi 
moves to the cell occupied previously by particle Pi+i. This behavior may be 
unwanted if the rule is supposed to express the motion of indistinguishable 
elements. The following motion representation has a 3-cycle; less trivial local 
cycles can be found in other PA (and it is even possible to construct PA with 
local cycles of arbitrary length). 

M12 = { OillO , OlflO , omo } 

Order preservation: We say that a PA preserves the order if there is no 
configuration in which a particle moves from position i to position i x , while 
another moves from j to j\i with i < j < j\ < h, or i x < ji < j < i . 
This behavior may be wanted, for instance, when modeling cars moving on a 
one-lane road. Order-preserving PA do not admit local cycles. 

Anticipation: We say that a PA with states {0, ... , q — 1} exhibits anticipa- 
tion, if for some configuration a cell in a state s receives, from other cells, a 
number of particles t such that t + s > q — 1 (as if, as a result of the rule, the 
neighboring cells "knew" that some particles will leave the cell). For q = 2, 
anticipation is what happens when a particle moves to a currently occupied 
cell. In the case of highway traffic, for instance, most drivers usually move 
anticipating the motion of the cars ahead (assuming that they are not going 
to stop), and this is considered in the models. An important case in which 
anticipation may be unwanted is the situation in which movement is optional: 
the particles may decide to stay in their place or to move, but if they move, 
they must obey the PA rule. In this situation, an anticipating rule would cause 
collisions (or overflows, if q > 2). 

Global cycles: We say that there is a global cycle if there is a chain of 
particles {pi}^, located at different positions, such that, for all i, particle pi 
moves to the cell previously occupied by particle Pi+\. The reason to call this a 
"cycle" is the following. By removing local cycles, the chain may be assumed to 
approach oo (or — oo) as m — > oo. If the configuration is — spatially — periodic 
of period p, then we may identify it with the torus Z p , and the chain is in 
fact a cycle. If the configuration is not periodic, as we follow the chain we 
will at some point find a repetition (due to the finite number of possibilities) 
of a block of length larger than \N\. At this point we may cut the part of 
the configuration starting with the block and ending with it, and repeat it to 
produce a periodic configuration with a cyclic replacement. A trivial example 
of global cycles is the PA that just shifts the configuration, as motion rule M 3 
given above. 

A special case are global cycles of anticipatory motion (which is always the 
case for global cycles when q = 2). Such cycles may be unwanted if we are 
modeling agents with local information, and we do not want them to use the 
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information "I am on a torus" , or "I am in an infinite queue" . 

Theorem 7 For any conservative CA F there exists a G E VA° such that 
n(G) = F and which preserves the order (and hence, has no local cycles). On 
the other hand, there do exist rules in CA° for which anticipation and global 
cycles are intrinsic: with any offset, they are not the projection of any PA 
without these features. 

Proof. The first part follows from the construction in Theorem 5: the canon- 
ical PA preserves the order. For the second, we just need to exhibit a CA for 
which the claimed property is true. 

Consider a CA F with the local rule / E £4° (2, 5) of Example 3. We will 
show now that any G E VA° such that 11(G) = F must have global cycles 
(and, in particular, anticipation), for any offset. Let G E VA°(q,n) be such 
that n(G) = F. Consider a configuration 

... 00000.111111... 1111111,0000000... 

where the sequence of l's is longer than 2n; the dot and the comma are there 
for reference. The image of this configuration, assuming an offset 0, is 

... 000011.111111... 1110011,0000000... 

If the particles in the middle of the configuration are moving, then they are 
moving without seeing any 0's; hence, the particles in the configuration 1 
(where all states are 1) would also be moving, and that would be a global 
cycle. On the other hand, if the particles in the middle are not moving, then 
there are two particles that are moving somehow from one end of the region 
of l's to the other, which is a contradiction, since the region is larger than the 
neighborhood of the PA. 

For any other choice of the offset, we obtain the same situation, except for an 
offset of 2. But in that case, we can consider the configuration 

000000.1010101010... 010101, 00000 

whose image is 

000000.0010101010... 010101, 01000 

and produces the same result as above: for any PA, there is an anticipation 
to the right, and it allows a global cycle. □ 

If anticipation and global anticipatory cycles may be intrinsic to a rule in 
CA° , then it is natural to ask about the decidability of these properties: given 
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a rule / G CA , can we decide if, for some offset, there does exist a PA without 
anticipation (and/or without global anticipatory cycles) from which the CA is 
the projection? The answer is not trivial. It is easy to check these properties 
on a given PA; however, to check them on a given CA rule, we must consider 
all the possible PA for which the CA is the projection. Consider the CA of 
Example 4: the canonical PA given by Theorem 5 shows anticipation, but a 
non-anticipating PA with the same projection does exist, and is described by 
the motion representation 

Mi 3 = { lib 010 } 

Hence, to find a non-anticipating PA we may need to drop the condition of 
order preservation. In fact, the situation is even worse. The canonical PA of 
Theorem 5 has a nice feature, stated in Corollary 6: a particle moves always to 
a cell that "sees" it in its CA neighborhood. As we shall see in the next sections, 
we can ask for the same feature when obtaining PA for monotone CA and for 
state-conserving CA; it seems to be "natural", and some authors have taken 
for granted that only motion rules with this property need to be considered. 
But to avoid anticipation, we may need to drop this condition too. It is violated 
in the PA described by M i3 , and it can be shown that anticipation cannot be 
avoided without violating it for the local rule of Example 4, combined with 
any offset. Thus, we must consider a rather large set of PA for a given CA 
rule in order to decide if anticipation can be avoided or not. 

Theorem 8 For a given f G CA° , it may be decided whether an offset £ exists 
for which the CA defined by f and £ is the 'projection of some PA without 
anticipation. If the answer is positive, that PA can be found. 

Proof. First we must notice that for a given / G CA° , there is at most 
one offset £ for which a non- anticipating PA can exist. This was already seen 
in the proof of Theorem 7 for the rule of Example 3: when we evaluated a 
configuration of the form 01 m 0, with an arbitrarily large m, an offset had to be 
imposed to prevent the movement of particles from one extremity of the l's to 
the other. For the general case, let us suppose that there is a non-anticipating 
PA G such that its projection is / with offset £ (and hence neighborhood 
{— £, ... , r}, with r = n—l—£). Then we may assume that G has a neighborhood 
{— L, ... , R} with L > £ and R> r. Consider a configuration c = 0(q — 1) M 0, 
with M arbitrarily large. We have the situation depicted in Table 3. 





A 


B 


C 


D 


E 


F 


G 


c 







(q-l) L - e 


(q _ l}M-(L+R) 


(q-l) R - r 


(q-l) r e 





c' 





f(!^ +r (q-l) t+r ) 


(g-i) L - e 


(q _ i)M-(i+fl) 


(q-l) R - r 


f((q - lY+rQl+r) 






Table 3 

Iteration on 0(q - 1) M 
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The PA must keep the particles in region D fixed: they only see cells in state 
q — 1 around them. If they move, then we would have anticipation (and a 
global anticipatory cycle) for the configuration where all cells are in state 
q — 1. Furthermore, the particles in regions B and C cannot move to region 
D (this would imply anticipation), nor to regions E, F and G (since M is 
arbitrarily large). Hence, G must match the L(q — 1) particles from regions B 
and C with the particles that d has in regions B and C, and this requires 

L(q-1) = (L-£)(q-l) + J2f(0_^0(q-l)...(q-l)) 

i=i - ' 

1 n—t 

which implies 

i n— 1 

/(Q Mj^i) ) (5) 

1 n—i 

Thus an offset is imposed for the existence of a non-anticipating PA. If the 
right side of (5) is not an integer, then we are done with /: any PA will have 
to show anticipation, and, moreover, a global anticipatory cycle (since the 
anticipation takes place through the arbitrarily large region D, and hence, 
will also take place in the configuration with all cells in state q — 1). 

We will attempt to define a non- anticipating PA G with neighborhood {— L, R}, 
where L = L\ + £, R = R\ + r, and Li, R\ are "large" numbers that will be 
precised below. If the attempt fails, we will show that a global anticipatory 
cycles must exist. 

The construction of G is inspired by the proof of Theorem 5. For every se- 
quence w = W-l, •■■ , Wr G Q l+r+1 , we will consider the configuration c = OwO 
with wq at the origin and will use it to define the new positions of the wo 
particles. As before, we write d = F(c), where F is the CA defined by / and 
the (now fixed) offset £. 

If G is to be non-anticipating, then there are certain cells were we already now 
that some particles cannot move. Without anticipation, the cell at position % 
can receive at most q — 1 — q particles from other cells. If c- > q — 1 — q, 
then at least d i — (q — 1 — q) particles in position i must stay there, without 
moving. For every i 6 Z, we define 

hi = (c- + Ci - q + 1)+ a» = Ci - h d, t = c ■ - h { 

Thus, (6j)i e z are the numbers of particles fixed at the different positions, (ai) ie z 
contains the numbers of particles in the preimage which need to be associated 
to some particles in the image, and (dj)jez has the numbers of particles in the 
image which need to be associated with some in the preimage. Clearly, a, = 
bi = ^ = Ofori ^ {-L l -£,R 1 +r}, and d { = for i £ {-L 1 -£-r, R 1 +r+£}. 
As in Theorem 5, we associate the two sets of particles, (aj)jgz and (dj)jgz, 



16 



from left to right. The PA is thus defined: the particles at the origin have 
found their new locations, possibly fixing some of them (bo), and associating 
the rest to the available particles in the image. Notice that in general this PA 
does not preserve the order of the particles. 

Suppose now that (for at least some w), a particle from the origin is sent to a 
position outside of {—Li,R{\. We consider the case in which it is associated 
to a position to the left of — Li (the other case is symmetric); this means that 
the non-fixed particles to the left of in c were not enough to match all the 
non-fixed particles to the left of —L\ in c' , i.e. 



-1 -Li-l 

Y a i < E d i 

i=-L x -i -Li-i-r 



-1 

E a i 


< 


-Li-1 -Li-1 

E d i- E a i 






i=—Li 




-L-L-l-r i=-L!-t 






-1 




-Li-1 -Li-1 


-Li-1 


-Li-1 


E a i 


< 


E E 


- E 6i 


- E a 

i=-L\-t 


i=—L\ 




-L-i-l-r i=-Li-£ 


i=-Li-^ 


-1 




-Li-1 -Li-1 






E a i 

i=—Li 


< 


E E ^ 

-L-i-l-r i=-Li-e 







(6) 



On the other hand, 

-Li-1 -Li-l+r 

E < < E < E (7) 

-Li-i-r -Li-t-r -Li-l-r 

(otherwise, the configuration Oc-^-e, ■■■ , C-^-i+r® would contradict / G CA° ). 
Combining (6) and (7) we obtain Y^T=-l-i a i < r ( ( l ~ l)i thus, by choosing 
L\ > r(q — 1)L 2 (where L 2 has still to be chosen) we can guarantee that there 
will be an interval I = {s, ... , s + L 2 — 1} C {—L\, ... , —1}, such that a* = 
for all % e J, i.e., all particles in this region are fixed. 

Now, we define t = 2max{£, r}, and choose L 2 = (t + l)g*. In this way, we 
guarantee the existence of a word u of length t that occurs at least twice, 
without overlap, in c s , ... ,c s+ l 2 -i- We can write it as u — u\u 2 , with = 
\u 2 \ = max{!,r}. Let — T and — T' be the positions to the left of the first 
occurrence of u, and to the right of its second occurrence, respectively. If we 
write x = c_ Ll _ £ , ... , c_ T , x' = c'_ Ll _ e _ r , ... ,c'_ T , y = c_ T /,^. ,c Rl+r , and 
y' = d_ T i,... ,c' Ri+r+l , we can rewrite c as c = OxuiU 2 u 3 u 1 u 2 yO for some u 3 , 
and we have that 

F(c) = F(0xuiu 2 U3Uiu 2 y0) = Ox'u[u 2 u' 3 u[u' 2 y'0 



17 



If we define c M = 0xui(u2UsUi) M i^yO, for an arbitrary positive integer M, we 
obtain that 

F(c M ) = F(0^ 1 (w 2 w 3 w 1 ) M w 2 yO) = 0xfu[(u' 2 u' 3 u' 1 ) M u' 2 y'0 

The situation is similar to the first paragraphs of this proof. Since M is arbi- 
trarily large, we have an arbitrarily large region were particles cannot be moved 
by any non-anticipating PA. If we choose M such that \uu 3 \M > L' + R! + 1, 
where {—I/, ... , R'} is the neighborhood of any candidate non-anticipating PA, 
we see that the PA will have to match the particles of xu\ with those of x'u[. 
This requires an equal number of particles in both of them, i.e., 

-T+t -T+t 

E c * = E < 

i=—L\—£ i=—L\—£—r 



-Li-l 


-T+t -T+f 


-Li-1 


-T+t 


E °i 


+ E a * + E & * = 


E 4 


+ E cj 


i=-L x -i 


i=—L\ i=—Li 


i=-Li-^-r 


i=— Li 


-T+t 


-T+t -Li-1 




-T+t 


E &*- 


E 4 = E 4 


- E c * 


- E a * 


i=— Li 


i=—L\ i=—L\—l—r 


i=-Li-< 


i=—L\ 


-T+t 


-T+t -Li-1 




-1 


E &*- 


E cj > E 4 


- E ^ 


- E °i > 


i=— Li 


i=— Li i=—L\—£—r 




i=— Li 



where the last inequality uses (6). We have arrived to a contradiction, since 
bi < €[, for all i. 

From all the previous discussion, we see that our construction of G, with 
neighborhood {-L x -l, ... , i?i+r}, L x = r{q-l){M)q\ R x = £(q-l)(t+l)q t , 
will move the particles from a position % to {% — L 1; ... , % + if not, then 
/ does not admit a non-anticipating PA, and it exhibits global anticipatory 
cycles, with any offset. 

We still have to show that G is really non-anticipating, and that F is its 
projection. For the first fact, we have to notice that the cell at position % is 
receiving at most di particles (it may receive less than di, since the assignment 
of non- fixed particles may fix some of them). For G to be non-anticipating we 
need 

di < q - 1 - Ci c ■ - bi < q - 1 - c { c • + q - q + 1 < 6, 

which follows from the definition of bi. 

To show that F is a projection of G, we follow again the scheme of Theorem 5, 
and consider the effect of adding an arbitrary word u to the left of w. Instead 
of c = 0u>0, we take now c = OuwO, and proceed as before: we define bi, di, di, 
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fix hi particles at position i, and associate the "free" particles from left to right. 
We will show that the destination of the Wq particles at the origin is exactly 
the same as before; this implies that applying G to a finite configuration c 
has the same effect of matching the particles with those from F(c), and we 
conclude that 11(G) = F. 

In order to show that the destination of the particles at the origin has not 
changed, we will show that the number of "free particles" sent from {— L ± , ... , — 1} 
to positions to the left of —L x is not changed. Since for % >— —L x we have not 
only Cj = Q, but also c[ = c- (and hence 6, = b iy di = a iy di = di), this implies 
that the destinations of the particles at the origin remain the same. 

The number of "free particles" sent from positions {—L 1 , ... , —1} to positions 
to the left of — L x is equal to the difference between the number of particles 
to the left of — L x in the image, and the number of particles to the left of — L x 
in the preimage. Thus, we have to show that 



-Li-l -Li-1 -Li-1 -Li-1 

E c 'i- £ Ci = £ d i- £ 

i=-L 1 -(-r i=-L!-e i=-L 1 -£-\u\-r i=-Lx-i-\u\ 

-Li-1 -Li-1 

E <= E 3-E«* (8) 

i=-L 1 -e-r i=-L!-e-\u\-r i=0 

Since F is conservative, and since = for i >= — L ± , we have 

Ri+r -Li-1 Ri+r+l 

E E E < (Q) 

i=—L\—i i=—Li—l—r i=—L\ 

and 

\u\-l i?!+r -Li-1 i?i+r+« -Li-1 Ri+r+l 

J2u t+ J2 ^= E 3+ E ^= E E < 

i=0 i=—L 1 -£ i=—Li—l—\u\—r i=—L\ i=—L 1 -£~\u\~r i=—L\ 

(10) 

We obtain (8) from (9) and (10). □ 

Corollary 9 The intrinsic presence of anticipation and the intrinsic pres- 
ence of global anticipatory cycles are equivalent properties for conservative 
CA rules. 



Proof. One direction is trivial, since global anticipatory cycles include an- 
ticipation. The other follows from Theorem 8: if the procedure fails to give a 
non-anticipating PA, then a global anticipatory cycle occurs. □ 
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A comment about the feasibility of finding a non-anticipating PA is due. Even 
if Theorem 8 asks, in principle, for the evaluation of a huge number of configu- 
rations ((? — 2)(g — l)( n-1 )te -1 )(* +1 )9 + n ~ 1 ) with t = 2max{f, r}), more practical 
implementations are possible: if we apply the procedure of the theorem to a 
configuration 0w- a , ... ,wo, ... , wpO, with a > £ and (3 > r, and we find that 
the wq particles at the origin move to positions in {— a + £, ... , /3 — r}, then the 
argument contained in the last part of the proof holds, and we do not need to 
evaluate the configurations of the form uw'v: the movement of these w will 
be the same in all cases. 

4-3 State- conserving CA 

Some CA satisfy a condition stronger that number-conservation: they conserve 
the number of cells in each state, throughout the iterations. We say that / G 
CA(q, n) is state- conserving if it verifies 

\w\-i H-i 

£ S a (f(w/w) k ) = £ 5 a (w k ) \JweQ*,aeQ (11) 

k=0 k=0 

where S a (x) is 1 for x = a and otherwise. It is plain to see that state- 
conserving CA are conservative, and that all elements in CA°(2,n) are state- 
conserving. State-conserving CA may be characterized as follows: 

Proposition 10 A necessary and sufficient condition for f e CA(q, n) to be 
state- conserving is that, for all a,x±, ... , x n e Q, 

S a (f(x 1: ... ,x n )) 

n-l 

= S a (xi) + <W(JV^0, X 2 , , X k+1 )) - 6 a (f (0^0, X U ... , X k )) 
^ ^ n—k n—k 



Proof. For each a G Q, we apply Theorem 2 to the density function 5 a . □ 

Notice that 5 a is a non-linear function. Thus, unlike the characterization of 
CA° in 3, this characterization cannot be used to create a linear system whose 
solutions would give all the CA satisfying the condition; however, it may be 
used as a test to look for state conservation in CA° . 

State conservation is local: If there is a state a in a configuration, then 
there must be a state a close to it in the image. Suppose this is not true. 
Then there are configurations where no a appears in the image of a window 
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of length L around the original a, for arbitrarily large L. We may choose L 
large enough to assure that there is a word of length n which is repeated in 
the window; we cut the configuration between the repetitions (keeping one 
of them), and obtain a periodic configuration without a state a in its image, 
which is a contradiction. 

The right PA for a state-conserving CA: If state-conservation is local, 
then the most reasonable way to look at the rule in terms of particles is to 
consider each state as a single particle, with different states corresponding to 
different types of particles (and, in fact, the numerical values of the states 
turn out to be irrelevant, and could be replaced by colors or letters). But 
then we would like to have a kind of particle automaton that takes a particle 
of type a and a surrounding configuration, and moves it to the position of a 
particle of type a in the image; it would be defined by a tuple (q, {—(-, ... ,r},g), 
with g : Q e x Q r — > N, and g would determine the new position of the particle 
currently located at the origin in N. If we want to keep the image of the states 
as occupancy numbers, then what we would want is a particle automaton that 
moves all the particles of a cell together: instead of evaluating to an arbitrary 
vector in N l , each g, L (in the notation for PA) would evaluate to a vector 

{hh - for so me 3 e N - 

V v ' 

i 

Constructing the right PA: Unfortunately, this is not the PA that will 
result from the application of Theorem 5 (unless we have q = 2, or a particular 
case like the shift). However, the construction given in the proof of Theorem 
5 may be fixed for the case of state-conserving CA: for each configuration, we 
assign now to each particle (i.e., to each state in each cell) the position of 
its correlative particle (from left to right) in the image of the configuration 
(this can be done, thanks to the state-conservation). The proof then proceeds 
exactly in the same way, and the resulting PA will have the same neighborhood 
(and the particles will arrive in the same zone inside it) as stated in Corollary 
6. In general, order will not be preserved. 

Example 5: Consider / G CA(S, 3) with code #6768185473053, and offset 
1. Using Theorem 5 we obtain the motion representation 

i 

M 14 = {21} 

However, this CA is state-conserving: a '2' will travel to the right as long as it 
is immersed in a background of l's. Depending on the application, it may be 
therefore more appropriate to use the special version of the construction (as 
described above), and obtain the motion representation 

2 

M 15 = {21,21} 
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where the numbers in the arrows might be dropped, since they will always 
represent the motion of the complete "particle", 1 or 2. □ 

4-4 Momentum Conservation 

So far we have considered one additive conserved quantity, the mass. It is 
natural to ask about other quantities that frequently follow conservation laws, 
as, for instance, momentum. Notice that this question does not apply to CA, 
but is natural for a PA. 

In a PA G = (q, N, (<?i)j = i v .. with N = {—£, ... , r}, we define the velocity 
of a particle at a given time step as the difference between its position at 
the next time step and its current position; equivalently, as the value that g.i 
assigns to it. Thus, the sum of the velocities of the particles at a cell i of a 
configuration c G S z is 

Ci 

V i( c ) = ^2[9ci(Ci- e , ... ,Cj_i,Cj+i,... ,Ci +r )] k 
k=l 

We will say that a PA preserves the momentum if, and only if, 

p(c)-l p(c)-i 
E ^( c ) = E W) VcGCp, (12) 

i=0 i=0 

where Cp are the periodic configurations of Q z and d = G(c). In other words, 
the function : Q e+r+1 — > Z defined by 

<f>(x ,... ,x e+r ) = Et^i^o,- ,x e -i,xe+i,... ,x e+r )] k 

k=l 

is asked to be the density of an additive preserved quantity for 11(G). Direct 
application of Theorem 2 yields the following proposition. 

Proposition 11 Momentum preservation is a decidable property of PA. 

In fact, in spite of being defined in terms of the particle representation, momen- 
tum preservation depends only on the conservative CA we are representing, as 
shown by the next theorem. Notice that this is not true for non-conservative 
CA. 

Theorem 12 Let F be a conservative CA and let G and G' be PA such that 
n(G) = II(G") = F. Then the following three are equivalent: 

(i) G preserves momentum 

(ii) G' preserves momentum 
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(Hi) F verifies 



EE c r c j = EE c 'j - c j Vc e °f, where c' = F(c), c" = F 2 {c) 

i£Z j<i i<=Zj<i 



Proof. Since momentum is an additive quantity, Theorem 1 implies that con- 
dition (12) can be tested on the configurations of Cp instead of Cp. Therefore, 
(i)<^(ii) follows from (i)<^(iii) together with the fact that 11(G) = IT(G"), and 
we just need to show (i)<^(iii). 

If F = U(G), then F(c) = G{c) for all c e Q z , and we can forget F. Let c be 
a finite configuration, c' = G(c) and c" = G(c'). Then what we must show is 
that 

EW = EW) EE<*-<4 = EE4--<? 

iez ieZ i& j<i ieZ j<i 

In fact, what we have is that 

Ci 

EE c r c i = E^( c ) = EEtfc^.- ,Ci +r )]j (13) 

ieZj<i i£l ieZj=l 

To see this, consider the contribution of each particle in c to each of the 
sums. On the left side, the j-th particle at contributes its displacement, 
v = [g Ck (c k _ e , ... ,Cfc_i,c fc+ i,... ,c fc+r )]j. Without loss of generality, suppose 
v > 0. The left side of (13) is the addition, over i E Z, of Y.j<iCj — c'j, the 
accumulated difference between c and c'. The particle is moving from k to 
k + v; thus, it contributes +1 to this sum, for the v terms corresponding to 
i — k, ... , k + v — 1, i.e., its movement contributes with v to the total sum. □ 

Corollary 13 Momentum preservation is a decidable property in CA° . 

Example 6: Most of the momentum preserving NCPA with small neighbor- 
hoods are trivial (identity, shifts); the rest of the cases consist of rules that only 
allow movements with zero sum, as in the following motion representations. 

Mi 6 = { 00110 , 00110 } , Mi 7 = {120} 

Of course, more sophisticated CA with momentum preservation can be con- 
structed for larger neighborhoods. □ 



5 The Monotone Case 



In this section we deal with the characterization and particle representation of 
monotone one-dimensional CA. In fact, we will talk almost exclusively about 
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non- increasing CA (CA~ ), but it must be noticed that this is equivalent to 
talking about non-decreasing CA (or, at least, each result about the former 
translates into a result about the latter). In fact, there is a one-to-one cor- 
respondence between the elements of CA~ and those of CA + , by replacing 
"particles" with "non-particles" and vice- versa: for each / G CA~ we have a 
/ G CA + defined by 

f(x 1: ...,x n ) = q-l - f(q - 1 - x 1: ... , q - 1 - x n ) 



5.1 An Only Sufficient Condition 



A first idea for a characterization of monotone CA would be the replacement 
of the equality for an inequality in condition (3), i.e., we would like to say 
that / is in CA~(q, n) if and only if, for all (xi, ... ,x n ) G Q n , 



n-1 



f(xi, ... ,x n ) < X! + J2 ' 



k=i 



f(0, ... , 0, x 2 , ... , x fc+ i) - /(0, ... , 0, X!, ... , x k ) 



n—k 



n—k 



(14) 



However, this condition is only sufficient. To see that it is sufficient, consider 
any periodic configuration, and add both sides of the inequality along a whole 
period: all the terms in the sums of the right side will cancel, and we obtain (1). 
On the other hand, the condition is not necessary, as shown by the following 
example. 

Example 7: Let F be the elementary CA 72. Here q — 2, the offset is 1, and 



f(x 1 ,x 2 ,x 3 ) 



1 if x 2 — 1 and x\ + x 3 = 1 
otherwise 



It is easy to see that / G CA" : the only way in which a cell can get a 1 in the 
image is to already have one in the preimage. On the other hand, condition 
(14) is not verified: 

i = /(i,i,o) > i + /(o,i,o) + /(o,o,i)-/(o,i,i)-/(o,o,i) = o 



□ 
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5.2 Domination and maximal elements in CA 



We will say that h G CA(q,n) dominates f G CA(q, n) if f(w) < h(w) for all 
w G Q n . It is easy to see that any rule / G CA which is dominated by a rule 
h G CA° belongs to CA' . This arises a natural question: are all elements of 
CA~(q, n) dominated by elements of CA (q,n)l In fact, this is true for q = 2, 
n = 3 (the elementary CA). However, the general answer is negative: 

Example 8: Consider / G £4(3, 2) defined by 



2 if xi = 2 

1 if X\ G {0, 1} and rr 2 = 1 
otherwise 



It is easy to see that / G CA" , since f(x\,X2) < xi. Now suppose that there 
is h G £4°(3,2) such that f < h. Since /i is conservative, it must verify 
h(0, 0) = 0, /i(l, 1) = 1 and h(2,2) = 2; since it dominates /, it must also 
verify h(2, 0) = 2 and h(2, 1) = 2. From (3) we have 



2 = h(2, 0) = 2 + h(0, 0) - /i(0, 2) , and 
2 = /i(2,l) =2 + /i(0, l)-/i(0,2) 



Since /i(0, 0) = 0, we get h(0, 2) = and hence h(0, 1) = 0, which is less than 

/(o,i) = i. 

The CA defined by / (with offset 0) keeps the 2's untouched, while the l's 
travel to the left until they hit a 2, and disappear (it is the projection of the 
PA described by M 5 in Section 3). Most small examples of non-dominated 
rules in CA~ follow the same pattern: a wall of some kind, and particles that 
move until they meet it and disappear. Intuitively, a CA dominating them 
would have to preserve everything which is being destroyed (and it cannot do 
this only "next to the wall"), at the same time that it follows the particles in 
their movement; therefore, it would have to increase the total mass. □ 

We will say that / G CA~ (q, n) is maximal if it is not dominated by another 
element of CA~(q,n). 



5.3 A Characterization of CA 



In the following proposition we show that monotony is decidable. Unfortu- 
nately, our characterization is computationally useful only for small values of 
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q and n. 

Theorem 14 Let f be a local rule in CA(q,n). Then f G CA~(q, n) if and 
only if 

\w\— 1 \ w \~ 1 

f(w/w) k < w ^ VweL(q,n) (15) 

fc=0 fc=0 

where L(q,n) is the set of all words in Q* such that 

(y^i mo d | w |, ... , W(i+ n —2) mod \w\*j 7^ (^j mod \w\i ■■■ j ^(j'+ra— 2) mod \w |) fori^j 



In other words, L(q,n) are the words that do not repeat a subword of length 
n — 1, when considered as circular words. They correspond to all the cycles 
in the de Bruijn graph B(n — 2, q) that do not repeat edges; their maximum 
length is that of the Eulerian paths in B(n — 2,q), which is g n_1 . 

Proof. Since (15) is the restriction of (4) to some configurations, the condition 
is obviously necessary. To show its sufficiency, we have to consider a word 
w G S* \ L(q, n). Then w must have a subword a of length n — \ which occurs 
twice in w. There are two cases: the occurrences of a overlap, or they do not. 

Case A: No overlap. Without loss of generality, we can assume that w 
begins with a, and write it as w = auav. In addition, we can assume that 
both era G L(q,n) and av G L(q,n): if not, we apply the whole argument to 
them (recursively). Thus, they verify (15), and we have 



\w\-l \w\-l 

f(w/w) k = (auav/auav) k 

k=0 k=0 

\au\— 1 |at)| — 1 

= Y (au/a) k + Y i av / a )k 

k=0 k=0 
\au\— 1 |c«j| — 1 | ^jj | — 1 

fc=0 fc=0 k=0 



Case B: With overlap. In this case, two occurrences of a overlap. It follows 
that a has a prefix (3 such that = fa mo d Note that 

f(w/w) = f{w/a) = f((3au/a) = f((3/a)f(au/a) 

As before, we can assume that inequality (15) is satisfied for the words au 
and (5 (in fact, since < n — 1, we apply it to f3 m , with m\/3\ > n — 1, and 
divide by m to obtain it for ft). Thus we get 
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M-i 

E f(w/w) k 

k=0 



|/3|-1 |cm|-l 

E (P/ a )k + E ( au / a )k 
k=0 k=0 
l/3|-l |cm|-l 

< E fa + E ( au )k 

k=0 k=0 



\w\-l 
E W k 



k=0 



□ 

As we said before, this test is not very practical, since it involves checking the 
condition for a large number of words, which grows as q q " . A careful listing of 
the cycles in de Bruijn graph B{n — 2, q) could make this a bit lower, but not 
much: there are at least (q\) qn ' g _(n_2 ) Eulerian cycles in B(n — 2,q). 



5.4 Particle Representation for CA 



Once we have characterized one-dimensional monotone CA, and since we know 
that one-dimensional conservative CA can be represented through particle 
automata, a natural problem to consider is the representation of non-increasing 
CA in terms of the movement of particles. Clearly, the projection of any PA 
is a non-increasing CA. For the converse, we will show how to construct a PA 
for a given non-increasing CA. 

More precisely, we will show how to associate to any cell in the image of a 
configuration, the location that its particles had in that configuration (this is 
a particular case of the particle identifications defined in [Kur03]). Consider a 
word w G Q n , with image f(w) > 0. Since the CA is non-increasing, we have 
that f(w) < YhZq u>i- We will impose that the f(w) particles in the image 
cell must come from the cells in which wo, ... ,w n -i were located. Moreover, 
the particles will be assumed to be contiguous in the preimage. Then their 
location in w is defined by a single value s(w), s(w) G {0, ... , X^o 1 w i ~ f( w )}i 
as shown in the following scheme: 



w 



(16) 



s(w) f(w) 



Denote by p(w,i) the number of particles contributed by each Wi to the cell 
that holds f(w) in the image; this can be better understood in the scheme 
above, by seeing how many of the f(w) selected particles fall in i, but its 
precise value can be written as 



p{w,i) = \ mm {s(w) + f(w),'Y w j} ~ max {s(w),y^w ? } (17) 

j<i j<i 



27 



and we have f(w) = Yh=q }p(w,i). 

Theorem 15 Let f G CA(q,n). Then a necessary and sufficient condition for 
f to be in CA~ is that there exists s : Q n — > Z, with s(w) G {0, ... , J27=o w i ~ 
f{w)} for all w = w , ... ,w n -i and s(w) = for all w with f(w) = 0, such 
that p(w,i) defined by (17) verifies 

n-1 

p({w n -i-i, - , w 2n -2-i),i) < w n _i for all w , ... , w 2n -2 e <2 (18) 

j=0 



Proof. 

Sufficient condition. Suppose that there exists s with these properties, and 
consider any finite configuration c G Q z . Then 



n-1 

E ^ c i' - ' c i+«-i) = E E P(( c i' - ' c i+"-i)>0 

n-1 

= HHP^ C 3-ii- ,Cj-i+n-l),i) < J2 C 3 



Necessary condition. Suppose that / G CA . For each w = w , ... , u> n _i in 
<5 n define s(iy) = if f(w) = 0, and 

n-1 f — 1 n+|u| — 1 

s(w) = 



«=0 I «=0 «=0 



E^ + mm{E^- E /HO- 1 )^ (19) 

otherwise. It verifies 



n— 1 n— 1 n— 1 

o < *h < E^-E/^/o"" 1 )* < E^-/H 

i=0 i=0 i=0 



where the lower bound follows from / G CA , and the upper one is obtained 
by taking v as the empty word (or as "0") in the definition of s(w). 

It is important to notice that s can be determined with a finite computation: if 
v is a word of minimal length where the minimum is reached, then \v\ < q n ~ x . 
The situation is similar to the proof of Theorem 14; since we are not taking 
circular words, we need a word of length greater q n ~ l + n — 1 to assure the 
existence of a repeated subword a of length n — 1, and this length is reached 
by wv. Then wv can be written as xayaz, and we can define v' such that 
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wv' = xaz, which verifies 

n-l Kl-l n+\v'\-l 

E + E < - E /(w/o- 1 )* 

i=0 i=0 i=0 

|aj/|-l |aj/|-l 

= s(«;)+ E E f( a V/ a )i < S ( W ) 

i=0 i=0 

where the last inequality is implied by the non-decreasing property: if it is 
false, then the periodic configuration with periodic pattern ay contradicts the 
monotony. Here we assumed that the two occurrences of a do not overlap; if 
they do, the reasoning is similar, akin to that in case B of Theorem 14. In 
both cases, what we get is a violation of the minimality of \v\. 

We still have to show that s(w) verifies (18). This follows from the fact that 
different particles in the image take different particles as their preimages. 
Consider a configuration c G Q z and two positive coordinates in F(c); we 
may assume that the offset is (i.e., F(c)o = f(co, ... ,c n _i)), and that the 
two positions are and some k > 0. Consider a scheme like (16), with the 
X^o™ 1 c i particles of the positions 0, ... , Cfc+„_i aligned in a row, and number 
them from left to right. Then the f(c , ... , c n _i) particles landing at position 
are taken as the particles labeled with numbers 

s(c , ... , c n _i) + 1, ... , s(c , ... , c n _i) + f(co, ... , C„_i) 

while the f(ck, ... ,Ck+ n -i) particles landing at k are taken as the particles 
labeled with the numbers 

k-l fc-l 

Cfc + S(c k , ... , Ck+n-l) + I, - , E °k + S ( C fc> - ' C k+n-l) + /(Cfc, ... , C fe+n _i) 
i=0 i=0 

All we need to show is that they do not overlap, i.e., 

fc-i 

s(co, ... , c n _i) + /(c , ... , c n _i) < J^c fc + s(c k ,... ,C fc+n _i) 

i=0 

Consider any word v G Q™, and put the word v' = (c n , ... , c k +n-i, Vq, ■■■ , v\ v \-i) 
in (19) for the definition of s(c , ... , c n _i). We obtain that 

n-l Kl-l n+K|-l 

S(C , ... , C„_i) < E C * + E U i - E /(COJ- jCn-lj^V "" 1 )* 
i=0 i=0 i=0 

fc-l fc+n-1 |v|— 1 fc-l 

= E C *+ E C i + E ^ - E/( C *'- >Ci+n-l) 
i=0 i=k i=Q i=0 

|u|+ti-1 

- E /(Cfc,... ^fc+n-l^/O' 1-1 )^ 
i=0 
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Thus, any v G Q* verifies 

k+n-l \v\-l \v\+n-l 

Ci+^Vi- /(cjfc, ... ,c fe+ „_i,t;/O n_1 )i 

i=k i=0 i=0 

k-1 k-1 

> s(c , ... , C n _i) - ^ d + Z( C *' - ' °i+n-l) 

i=0 i=0 

and therefore, taking the minimum over all v, 

s(Cfc, ... , Ck+n-l) 

k-1 k-1 

> S(C , ... , C„_i) - ^ Cj + XI f( C ii - > C i+n-l) 

i=0 i=0 

fc-1 

> s(c , ... , c„_i) -^c« + /(c , ... , C„_i) 

i=0 

□ 

Corollary 16 Let F be a CA with local rule f G CA~(q, n). Then there exists 
G G VA(q, 2n - 1) such that 11(G) = F. 

Proof. The functions s and p as in Theorem 15 determine the particle au- 
tomaton. For any configuration wq, ... , w n -i, ... , W2 n -2 with w n -\ > 0, we define 
the motion of the u>„_i particles as follows: the number of particles moving 
from n - 1 to i, i = 0, ... , n - 1, is given by p» = p((w n _i_j, ... , w 2n -2-i),i), and 
the number of particles that die (move to f) is given by w n -i — J27=o Pi- m 
other words, we define a g Wn _ x such that 

: [g Wn -i ( W 0, - , Wn-2, W w ... , W 2n -2)}k = i} =Pi 

n-l 

and #{£; : [^.^Wo, ... ,w n - 2 ,w n , ... ,w 2n -2)]k = = - 

i=0 

□ 

Corollary 17 (Another test for CA monotony) The previous theorems 
establish the equivalence between (one- dimensional) non-increasing CA and 
particle automata. This provides an alternative characterization of CA~ : to 
see if a certain CA belongs to this class, check the existence of a function s 
with the properties of Theorem 15. Finding it in the way described in the proof 
would need (in definition 19) as much work as the characterization in Theorem 
14; it is usually easier to check for its existence by testing the possible functions 
s. Since s : Q n — > Z, and < s(w) < J2 w i ~ f( w ), an upper bound for the 
number of possible s is (qn) qn , which is more than the work in Theorem 14, 
but using that s(w) = 0, p(w,i) = for f(w) = 0, and applying the necessary 
conditions of Theorem 15, the possibilites for s can be drastically reduced in a 
practical implementation. 
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Remark: There is no canonical PA. As seen in Section 4, we can always 
take a "canonical" particle representation for a (one-dimensional) conservative 
CA, which is unique, and is characterized by the preservation of the order of 
the particles. In the monotone case, there is no canonical form: for instance, 
when two particles are close to each other and one disappears, there is an 
arbitrary decision favouring the survival of one of them. 
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2 
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89 


? 


? 



Table 4 

Some demographics: Number of one-dimensional conservative, state-conserving, 
non-increasing, and maximal non-increasing rules, for some q and n. 



6 Conservation Associated to Blocks 

In [Mor03], Theorem 6, a method was described to decide, for a given CA, 
whether its states can be relabelled with integer numbers so as to make the 
CA number-conserving; in Section 4, we see that this can be used, in turn, 
to see if the dynamics of the original CA can be understood in terms of an 
operator acting on a system of indestructible particles. 

In this section we consider a generalization of that theorem, giving values to 
the words of a given length, instead of the single states; we want this function 
of the blocks to be preserved by the iteration of the CA. In other words, we are 
looking for an additive conserved quantity; however, we will impose a further 
condition: it must separate the images of the blocks with different values at 
the origin. 

We want to view the function as an assignment of a value to the state in a 
cell, but a value which depends on the states of its neighbors: the same state 
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a may represent 2 or 3 particles, depending on whether or not it is followed 
by, say, another state a. Thus we can detect if the given CA may be seen as 
the projection, to fewer states, of a conservative CA which in turn is seen as 
an interaction of particles, restricted to certain configurations. In particular, 
this may automatically detect particle-like behavior in some non-conservative 
CA. 

Theorem 18 Let F be a one dimensional CA with state set ScZ, and let 
b be a positive integer. Then it may be decided whether or not there is a b- 
neighborhood- dependent relabelling of the states of S, which distinguishes the 
states, and whose sum is preserved by F. 

Proof. Consider a cellular automaton F with states ScZ, neighborhood 
size n, and local rule / G CA(q, n). Without loss of generality, the offset can 
be assumed to be 0. Let be a relabelling map, <fi : S b — > Z. For c G S z , we 
define $(c) as 3>(c)j = ... , c i+b -i). Let be the induced CA that makes 
the diagram 

Fi F^ I 
S % *> $(^) 

commute. F^ acts on the subshift Q(S Z ), and has the same neighborhood as 
F. We want to determine possible mappings (p that would make F^ number- 
conserving. In other words, we want to find (f> such that 

p(c)-l p(c)-l 

<f>(Ci, ... , C i+6 _i) = ^ 0(/(Cj... Ci+n-i),... ,/(Ci +6 _i... Q +6+n _ 2 )) (20) 
fe=0 fc=0 

for all c G Cp. Fix an arbitrary a G S. From Theorem 2 we have that a 
necessary and sufficient condition for (20) is that 

0/(a o , ... , a b+n ^ 2 ) - <j>(a , - , a 6 _i) = 

n+6-2 

{ - 0/(«> - > « 5 °o, - , a n +6-2-i) + (j>f(a, ... , a, ai, ... , a n+6 _i_j)} 



i=i 

6-1 



(21) 



+ Y <K a > - > a o, - , a b _j_i) - 0(a, ... , a, ai, ... , a fe _;) 



i=i 



holds for all a , ... , a n +6-2 G 5; here we use the notation 

(f>f(x Q , ... , x b+n - 2 ) = (f>(f(xo, - , x n -i), ... , f(x b -!, ... , x b+n _ 2 ))- 

This condition can be used to get an homogeneous linear system, where the \ S 
unknown values are {<p(w),w G S b }. The solution will be a linear subspace 



b 
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V C WL sb ; if V = {0}, then the algorithm gives a negative answer. If V ^ 
{0}, we are not ready yet: we demand the solution to distinguish between 
different states in the first position (we are interpreting the image of the block 
ao, ••• , as the value assigned to a ). Thus we ask that 

0({a} x S^ 1 ) n <p({(3} x S^ 1 ) = Va ^ p. 

or, in other words, cannot belong to the hyperplane (p(u) = 0(f), for u, v G 
S*, Mo 7^ fo- Let {-E u ,v}(n,i))G/ be the collection of these hyperplanes, where 
/ = {(u, v ) : u, v G S b , Uo 7^ t>o}. Then we want to see if V \ (U( u ,v)eiE U!V ) ^ 0. 
Since each V fl E u>v is a subspace of V, and linear spaces cannot be finite 
unions of proper subspaces, we have: 

v \ ( |J £ u>w ) = ^ v = (J vn£ v 

(u,u)e/ (u,v)ei 

3Me;:V = yn^„ « 3(u, u ) G / : 1/ c s tt)t , 

and this last condition can be checked by adding the hyperplane equations to 
the linear system. 

If solutions do exist, then there are integer solutions and they can be found 
in finite time. Any such solution will be the desired assignment of values to 
the states in the original CA, that depend on their contexts, and whose sum 
is preserved by the dynamics of the CA. □ 

The range of values for b. There is not always a value of b that allows us to 
find a number-conserving relabelling of the states: this is the case, for instance, 
of the trivial CA with states {0, 1}, neighborhood {0}, and /(0) = f{l) = 0. A 
natural question is: For a given CA, how many values of b must we check before 
we conclude that no b will work? This is an open question at the moment of 
this writing, since we do not currently have a bound for b; intuitively, we can 
expect that such a bound exists, since there is no reason why the considered 
neighborhood should be arbitrarily large with respect to the neighborhood of 
the CA. The problem is akin to that stated (and left open) first in [Tak89] and 
then in [HT91] of bounding the range of the additive conserved quantities for 
the elementary reversible CA considered there; both are particular cases of the 
general problem of bounding the range of the additive conserved quantities of 
a given CA. 

Particle representation for the relabelled CA. The construction given 
in Theorem 5 can be trivially restricted to any subset of the possible config- 
urations of states; in particular, it can be restricted to the subshift $(S' Z ), 
yielding a particle representation of the operation of our CA on the states and 
its neighbors. 

Example 9: Consider / G CA(3,A) with code 64056. It is not conservative, 
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Fig. 2. Evolution of the rule from Example 9, with colors associated to 0, 1 and 2 
particles. 

and it can be checked that there is no reassignment of values to its states that 
makes it number-conserving. However, there are solutions of the system 
given by equation (21), for b = 3. The smallest such (in norm) is 

0(0**) = 
0(100) = 2 
0(101) = 0(110) = 0(111) = 1 

We obtain for it the motion rule 

2 rv rx 

2000 2001 2002 10 111 112 

restricted to the configurations where no 1 is followed by two zeros, and all 2 
are. This is a representation of the original CA in terms of the interaction of 
indestructible particles, and highlights one of the propagating defects in the 
evolution of this CA. □ 



7 Conclusion 

This paper presented particle automata as simple systems of interacting par- 
ticles, that give a complementary view on the behavior of conservative and 
monotone cellular automata. They are defined in terms of a global opera- 
tor induced by a local rule with some finite neighborhood (in analogy to the 
CA formalism), and may be easily extended to higher dimensions or different 
topologies, as we expect to do in a future article. 

The special class of conservative PA is equivalent to the class of conservative 
CA; we discussed several properties that the conservative PA may exhibit, 
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some of which depend only on their CA projections. These intrinsic proper- 
ties of the particle representation include anticipation, the existence of global 
anticipatory cycles, and momentum conservation, in a natural sense defined 
here. Theorem 8 shows that the intrinsic presence of the first two is equiva- 
lent, and that this presence is a decidable property. Some properties that are 
not intrinsic are the preservation of order and the existence of local cycles, 
both of which can be always avoided by taking the "canonical" PA given by 
the equivalence result. Momentum conservation is an interesting case, since 
it turns out to depend only of the conservative CA, in spite of being defined 
(and making sense) in terms of a PA. 

The conservative PA given by the construction algorithm may fail to be the 
most useful one for a given CA. This is the case with state-conserving CA; for 
them, the natural particle representation is one where the different states rep- 
resent different particles (instead of occupancy numbers). We gave a character- 
ization of state-conserving CA, and explained the construction of appropriate 
PA for them. 

We have extended the work on conservative CA in two directions: on one 
hand, we have considered CA that are not conservative, but nearly so, since 
we can see them as the projection of particular configurations of a conservative 
CA with more states; their dynamics can be then understood in terms of 
indestructible particles. 

On the other hand, we have considered monotone CA, and we have shown that 
they can be characterized, and, moreover, that they can also be represented 
by means of a particle system. This equivalence is less obvious than in the 
conservative case; one of the difficulties is that no "canonical" PA exists for 
monotone CA. 

One open problem was already mentioned in the previous section: in Theorem 
18, to find a bound for the block size b that needs to be considered before 
the existence of a solution is discarded. Another natural problem, which 
arises from the text, is the following: to extend Theorem 18 to allow to be 
monotone. This would be very helpful, since it would help to detect particle- 
like behavior that includes decay or annihilation. Such an extension is already 
possible with the current results, but it would be computationally hard; the 
challenge is to make it in an efficient way. 

Please notice that some additional information, like rule lists, as well as some 
software, is available at our website: 

http : //www . dim . uchile . cl/ ~anmoreir/ncca/ 

Other software (like the C++ routines used to determine and solve the linear 
systems for listing the CA rules) is available, by request, from A. Moreira. 
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